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1 Historical Introduction 



Much phenomenology in QCD requires the use of the operator product expansion [1-3] and many 
generalizations such as 'factorization theorems' [4]. Among the ingredients are matrix elements of 
particular gauge-invariant operators, which correspond to parton densities (or distribution func- 
tions). The properties of these operators under renormalization are vital to all QCD calcula- 
tions, and one serious complication arises because gauge-invariant operators mix 1 with certain 
gauge- variant (non-gauge-invariant) operators. The renormalization directly determines the phe- 
nomenologically important anomalous dimensions of the operators — generally used in the form of 
Altarelli-Parisi splitting coefficients. 

The extra operators that mix with the gauge-invariant operators are unphysical — we will call 
them 'alien' operators. It has been known since the earliest days of QCD that one must demonstrate 
that these alien operators do not contribute to physics. Three theorems apply to the decoupling: 
One is that a basis can be chosen such that the alien operators are BRST-exact. Next, physical 
matrix elements of BRST-exact operators are zero. The last theorem is a trivial consequence of 
the second: The renormalization mixing matrix is triangular — alien operators do not mix with the 
physical operators. The theorems to establish this have been proven in their strongest form by 
Joglekar and Lee [5] and more recently by Henneaux [6]. 

Unfortunately, recent calculations by Hamberg and van Neerven [7, 8] contradict these general 
theorems. Their results therefore throw into doubt the basis of all higher order perturbative QCD 
calculations. Our purpose in this paper is to resolve this contradiction between the theorem and the 
calculations. We will show that the contradiction is only apparent, and that it arises from certain 
subtle infra-red (IR) problems that are unluckily intrinsic to the usual algorithms for doing pertur- 
bative QCD calculations. However, the problem of efficiently performing practical calculations is 
left for future work. 

The immediate motivation for the calculations by Hamberg and van Neerven was a long-standing 

1 Multiplicative renormalization is not sufficient to remove the infinities from Green functions of arbitrary composite 
operators; counterterms corresponding to different operators are needed. 
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discrepancy between calculations of the two-loop anomalous dimensions of the twist-2 covariant 
gluon operators in Feynman gauge [9-11] and the light-like axial gauge [12]. Since these anomalous 
dimensions are measurable, calculations performed in different gauges should agree, and this can 
readily be shown by the methods of [13], provided that one assumes the Joglekar-Lee theorem. 

Hamberg and van Neerven repeated the Feynman gauge calculation and discover that the older 
calculations [9-11] are in error because they assumed the applicability of the theorem that the 
renormalization matrix is triangular. Hamberg and van Neerven show that the renormalization 
matrix appears to be non-triangular. Their calculation supports the otherwise dubious light-cone 
gauge result and is in accord with supersymmetry arguments [8]. 

The roots of this failure are already present in the one-loop part of the calculation. Although 
Hamberg and van Neerven do not remark on it, their calculation shows that the finite part of 
a physical matrix element of their alien operator is nonzero at one-loop order, contradicting the 
second of the Joglekar-Lee results mentioned previously. They perform their calculation for the 
whole tower of twist-2 covariant gluon operators, but the problems are present for the simplest 
of these operators, the energy-momentum tensor 0^, for which the renormalization theory was 
worked out by Freedman, Muzinich and Weinberg [14,15]. The alien operators in that analysis are 
not manifestly the same as those used by Hamberg and van Neerven. The form of 6^ given by 
Freedman et al. is in agreement with the general theorems of Joglekar and Lee and of Henneaux. 
However, the gauge- variant operators used by Hamberg and van Neerven are obtained from the 
analysis of Dixon and Taylor; it is not evident that these operators are BRST-exact. 

This is where we start: A sufficiently detailed analysis of the energy-momentum tensor at one- 
loop order is enough to locate the source of the contradiction. We will verify that at one-loop 
order, the renormalization given by Freedman et al. is in fact correct. However, the one-loop 
gluonic matrix element of the alien operator fails to vanish at zero momentum transfer. 2 ' 3 We will 

2 The momentum transfer is denned to be the sum of the momenta flowing into the inserted operator vertex. 

3 Harris and Smith in a recent preprint [16] calculate a nonzero gauge dependence for the gluonic matrix element 
of the gauge-invariant part of at nonzero momentum transfer. This is again in contradiction with the general 
theorems. 
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find that the source of this incongruity is an infra-red divergence, but the divergence is not in the 
calculation of the matrix element. Rather, it is a quadratic divergence in the proof that the matrix 
element vanishes. 

The source of the divergence makes it clear that the proof of the theorem on the vanishing of 
the alien operators should be correct when one applies it to physical states. The problem arises 
when one considers matrix elements in an off-shell gluon state and then takes the gluon on-shell. 

But this clearly threatens the rationale for the usual methods of doing perturbative QCD 
calculations. Moreover, the renormalization matrix that Hamberg and van Neerven calculated and 
found to be non-triangular presumably includes some infra-red renormalization, contrary to what 
should be done. 

In Section 2, we state the Joglekar-Lee theorems. In Section 3, we list our conventions for 
pure-gauge Yang-Mills theory. In Section 4, we give the results of the one-loop calculation of 
two-point Green functions with the energy-momentum tensor operator inserted at zero momentum 
transfer and derive the renormalization mixing matrix. The calculation at nonzero momentum 
transfer is currently underway. The off-shell results, as well as the physical matrix elements, will 
be reported in the near future. The Appendix contains: a brief discussion of 'right derivatives'; 
the full Lorentz tensor structure of the two-gluon Green functions, which are abbreviated in the 
text, with a separation into leading- and higher-twist pieces; a list of the Feynman graphs with 
composite operator insertions used in the calculations; and the Feynman rules for the operator 
vertices considered in this paper. 

2 Renormalization of Gauge- Invariant Operators 

In this section, we state the three theorems that Joglekar and Lee [5] proved on the renormalization 
of gauge-invariant operators. In [13], the theorems are stated and all the easy parts are proven. 

Let G{ denote a set of gauge-invariant operators that mix under renormalization, and let A{ 
denote the set of alien operators with which they mix under renormalization. (We define 'alien' to 
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mean 'not gauge-invariant'.) Finally, let E{ denote the set of operators that vanish by use of the 
equations of motion 4 and with which the previous two sets of operators mix under renormalization. 

The first of the Joglekar-Lee results is that the basis of the alien operators A{ that mix with 
gauge-invariant operators can be chosen so that they are all BRST-exact, i.e., they can be written 



as 



2.1 



where we will call B{ the 'ancestor' of A{. 

The second theorem is that physical matrix elements of the BRST-exact alien operators, 
<Wst-B 8 , are zero. 

The last of the theorems is that the renormalization mixing matrix is triangular: 
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/ 
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Zaa 
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Of these theorems, the hardest to prove is the first. It can easily be shown that all counterterms 
to BRST-invariant operators are themselves BRST-invariant [13]. Then one must prove that any 
BRST-invariant operator is a linear combination of gauge-invariant operators and BRST-exact 
operators. Up to operators that vanish by the equations of motion, this is supposed to be proven 
by Joglekar and Lee [5], but we find that proof very hard to understand. A simpler proof on the 
basis of cohomology theory is presented by Henneaux in [6]. 

A simple proof of the last two theorems is given in [13]. The vanishing of physical matrix 
elements of the alien operators follows from a simple Ward identity involving the BRST variation 



Matrix elements of E t must vanish, but Green functions of E t do not. 
5 We use the convention of [17], where the wavy equal sign means that the relation is only true after one or more 
of the equations of motions have been used. 
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(also called a Slavnov- Taylor identity), once one knows that only BRST-exact operators are needed. 
This result trivially generalizes to show that Green functions of these alien operators with BRST- 
invariant operators are zero. BRST-invariant operators include gauge-invariant operators and the 
BRST-exact operators that comprise all our alien operators. 

The third theorem, on the triangularity of the renormalization matrix, immediately follows [13]. 
If any on-shell, physical matrix element of an unrenormalized operator in class A is to vanish, then 
its pole piece must also vanish on shell. 6 Since at least some of the physical matrix elements of 
any gauge-invariant operator are non- vanishing, it follows that the entry Zag must be zero; no 
operators in class G can mix with the operator in class A. Similarly, Zeg an( i ^EA m the bottom 
row of the mixing matrix must be zero because an unrenormalized operator in class E must vanish 
by the equations of motion, therefore its pole part must also vanish by the equations of motion. 

Note that at the level of pure Feynman graph calculations, a physical matrix element is one 
with the gluon polarizations being purely transverse and with the states being on-shell quarks or 
on-shell gluons. 

Prior to the work of Joglekar and Lee, it was shown by Freedman, Muzinich, and Weinberg 
[14,15] how to construct a finite energy-momentum tensor for gauge theories. Their operator can 
readily be seen to satisfy the first Joglekar-Lee theorem, as we will explain later. 

The problem we now face is that the calculations by Hamberg and van Neerven appear to violate 
all of the above theorems. 

3 Yang-Mills Conventions 

In this section, we list some common objects in Yang-Mills theory to exhibit our conventions 
and notation, but also because some, such as the energy-momentum tensor, play a pivotal role 
throughout this article. 

6 'Pole' in this context means a singularity as the dimension of space-time is varied. We are assuming the use of 
minimal subtraction with dimensional regularization to perform the renormalization (see Section 3.5). 
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3.1 Lagrangian Density 

The effective Lagrangian density of pure-gauge Yang- Miffs theory in generai covariant gauge is, in 
terms of unrenormafized (bare) fiefds and parameters (designated by hats), 

jC(x) = -lpr(x)F^ a (x) - f\[8 ■ A a (x)][d ■ A a (x)} + [^if^^f^fi)],, (3.1) 

where the antisymmetric fieid strength tensor is given by 

Pr(x) = d»A" a (x) - d"A£(x) - gc abc AZ(x)A" c (x) (3.2) 

and the covariant derivative acts on fiefds in the adjoint representation of the group as foffows 

[D ^{x)Lb{x)] a = Df, ac (x)u c (x) = [d^Sac - gCabcAf, b (x)]u c (x) = d fl u a {x) - gCabcAf, b (x)u c (x). (3.3) 

We are defining the Grassmann fiefd f) a (x) to be the anti-ghost and the Grassmann fiefd u a (x) to 
be the ghost. The parameter g is the coupiing strength, c abc are the structure constants of the 
underfying Lie aigebra SU(N), and A is the arbitrary gauge-fixing parameter in generai covariant 
gauge. The coior indices in the adjoint representation a, 6, c, . . . range from 1 to N 2 — 1. 

3.2 Euler-Lagrange Equations of Motion 

The Eufer-Lagrange equations of motion, using right derivatives for the Grassmann variabies 7 are 
where 

$ a G {A va ,u a ,fia}. (3.5) 



7 See Section A.l in the Appendix for a discussion of right derivatives. 
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We have 



(D^) a + \d v d ■ A a + gc abc (d v fj b )u c = 0, (3.6a) 
(D^f,) a = 0, (3.6b) 
d^b^) a = 0. (3.6c) 



3.3 BRST Symmetry 

The gauge-fixed effective Lagrangian density is not gauge-invariant but is invariant under the 
foiiowing giobai symmetry [18] 

S BRST Lj a = -lgc abc u b u c 8{, ( 3 - 7 ) 
<WsT??a= A<9 • A a 8£. 

Here, is a constant parameter with Grassmann parity 1, that is it anticommutes with the 
(anti)ghost fieid components (and the fermion fieid components, if there were any), but commutes 
with everything eise. We introduce the notation 

in anaiogy with the right derivatives for Grassmann variables to mean that is commuted or 
anticommuted to the extreme right and then removed. This variation, called the BRST variation, 
is a symmetry of the Lagrangian density, Eq. (3.1), since the change in the Lagrangian density is 
a four- divergence, without invoking the equations of motion 

S -k^jr = -Xd^(D^) a d.A a ]. (3.9) 
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The important property of nilpotence, 



filR ST ( An y thin 9) = o, (3.10) 

holds only after using one of the equations of motion, Eq. (3.6c), which will be called the 'trivial 
equation of motion' in what follows. 8 

3.4 Energy-Momentum Tensor 

The symmetric, conserved energy-momentum stress tensor density can be constructed from the 
canonical tensor by using Belinfante's procedure [20,21]. This is also the tensor proven by Freed- 
man, Muzinich, and Weinberg [14, 15] to have finite Green functions with renormalized external 
fields. It is 

—-a C - F F p 

-g^dp(A p a d ■ A a ) + A(i„ a d v d ■ A a ) + \{A V a d^d ■ A a ) ( 3 - n ) 

The gauge-invariant piece is 

a(GI) _ 1 p P 7T p _ p p p (3 121 

8 If the Lagrangian formulation with the Nakanishi-Lautrup field is used, as in [18] and [19], no equations of motion 
are needed to demonstrate the nilpotence of the BRST variation. 
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The gauge- variant piece is everything else, 





g^d p (A p a d ■ A a ) - \\{d ■ A a )(d ■ A a ) + (dPfi a )(Dpu;) a ] 



(3.13) 



+X(A lia d v d-A a ) + X(A va d ll d-A a ) 



The gauge- variant piece of the energy- momentum tensor is the BRST variation of an 'ancestor' 
operator 



and the last term vanishes by the trivial equation of motion, Eq. (3.6c). The ancestor operator 
defined above is not finite, that is it does not have finite Green functions with all renormalized 
fields. In Section 4.4 we present the finite ancestor operator after introducing the 'renormalized 
BRST variation'. 

The BRST variation of the GV piece vanishes (up to the trivial equation of motion) because 
of the nilpotence of the BRST transformation. We say that a gauge- variant operator such as 
0^t V ^ is BRST-exact if it has an ancestor. The BRST variation of the GI piece vanishes without 
using the equations of motion because the BRST variation of the gluon field is based on the gauge 
transformation of that same field, so any gauge-invariant quantity is automatically BRST-invariant. 



an^stor = {d v f, a )A^ a + {d^f, a )A v a - g^[\f, a d ■ A a + (d p fj a )A p a ], (3.14) 



since 




(3.15) 
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3.5 Renormalization 

We use multiplicative renormalization 



I 

A — Z 2 A 



l 



fi — 7 2 r> (3.16) 

A= Z\X (It is known that Z\ = Z^ 1 to all orders) 
9= Z g g\i e 

and dimensional regularization in 4 — 2e space-time dimensions with the Modified Minimal Sub- 
traction (MS) scheme. 9 See [13] for a thorough treatment of the subject. 

The Lagrangian density, Eq. (3.1), can be written in terms of renormalized fields and parameters. 



3 Hamberg and van Neerven work in 4 + e dimensions. 
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This is the same Lagrangian density so the same symbol, £, is used to represent both quantities 10 

£ =-\F^F^ v a - ±A(3 • A a )(8 ■ A a ) + (d^ a )(D^) a 
-\6Z n (d^A v a - 8 V A^ a ){d»A v a - d v A^) 
+ lSZ III gii e c abc (d l _ t A 1J a - d v A tl a )A^A v c 

-\SZ Iv g 2 ^c abc c ade A^ h A v c A$A v e ( 3 - 17 ) 
-\(Z x Z A -l)\(d-A a )(d-A a ) 

+6z (d^ a )d^ a 

-6Z I gn t c ahc (d^ri a )Alu c , 

where 

i i 

7 — 72 72 
^0 = 

III 

7 — 72 72 72 ry 

l = ^ A ^ ^ 9 

Z u = Z A (3-18) 

3 

z iu= z l Z g 

Zj iv= ^A^g 

Notice that the renormalized coupling g is dimensionless and that we have introduced a param- 
eter [i with the dimensions of mass. The renormalization constants 6Zi = Z{ — 1 have been given 
Roman subscripts which label the number of gauge fields in the counterterm vertex. 

Since the coupling is universal, the three different interaction vertices have associated renor- 



3 Some authors use the term 'renormalized Lagrangian density,' but it is not always clear what is meant. 
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malization constants related by the following 'renormalization constant Ward identities' 



Zg - Z I Z A 2 Z Q 



Z iii Z A 2 



3.19 



In the Minimal Subtraction (MS) renormalization scheme, the counterterms are the negative of 
the pole part only, with no finite component. In the Modified Minimal Subtraction (MS) scheme, 
the ubiquitous Euler's constant, 7^, and the natural logarithms of 4ir are absorbed into a new 
renormalization mass parameter, j2, defined by 



H = ji 



( 4vr 

\ e lE 



(3.20) 



Applied under dimensional regularization, the counterterms in either scheme are proportional to K 
We list the multiplicative renormalization constants in the MS scheme to the order needed for 
an 0(g 2 ) calculation of the Green functions considered in this paper 



Z = 1 



Z A = 1 



Zx= 1 



1 9 2 



167T 2 



1 9 2 



167T 2 



C A 



C a 



1 9 2 



e 167T 



Z g = 1 + 0(g*) 



where C(A) = N for the gauge group SU(N). 



0(9 A ) 
0(9 A ) 

0(9 A ) = Z~/ 



(3.21) 



3.6 (Modified) LSZ Reduction 

The residue of the propagator pole is used in the LSZ (for Lehmann, Symanzik, and Zimmermann) 
reduction formula to derive the S-matrix from Green functions. The basic idea is that the S-matrix 
is obtained from the asymptotic behavior of Green functions for large times (t ±00), and this 
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behavior is governed by the singularities of the external propagators. We use a modified version of 
this procedure to handle the infra-red divergent logarithms that appear in this massless theory. 11 
Let be the self-energy. Then 



gluon 



(P)8ab 




3.22a) 




(3.22b) 



where the cross-hatched blob represents all one-particle irreducible amputated graphs, including 
counterterms so that the blob has no ultra-violet (UV) divergences. We isolate the p 2 dependence, 
defining II (p 2 ) by 



KLon(p) = (p 2 9 aT -p a p T Mp 2 ), 



(3.23) 



noting that the gluon self-energy is purely transverse due to a Ward identity [13]. 

In a massless theory, the singularities in propagators, as a function of p 2 , are not simple poles 
after higher order corrections are included. 

The dressed propagators are then 



KKp) 



p 2 + ie 



1 



1 + n(p 2 



(7 T 

p p 
p 2 + ie 

> 



and 



V ah {p) 



n(p 2 



i + n( P 2 ) 

gluon9 aT ^a.bC 2 gluonP a P 7 



p 2 + ie 



(p 2 + ie) 2 



(3.24a) 



(3.24b) 



p 2 - Y, ghost (p) + ie p 2 + ie 

where c 2 is the residue of the propagator pole (the coefficient of p 2 in the denominator). To one-loop 



'We do not have a complete justification of our algorithm. 
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order, the gluon self-energy is 



gluon 



(p) 



167T 2 



Ca [P 9aT ~ PaPr) < 



4 V 1 A 



1 - 4 



-21 + Mn 

9 T 3 



1 + iln 



+ 0(g% 



(3.25) 



which implies that 



n(p 2 



167T 2 



CU4 



-21 + Mn 

9 T 3 



iln 



+ 0(g% 



(3.26) 



and the ghost self-energy to one-loop order is 



^ghost(p) 



167T 2 



Cap 2 



In 



P 2 



4 V a; \ 



<% 4 ) 



(3.27) 



The singularities in the propagators are not simple poles, but the leading power, with logarithmic 
corrections, is governed by the large-time behavior of the propagator. So to define the residue we 
use the following formulae, which would be valid when the physical masses are nonzero: 



gluon 



^{p 2 [i + n(p 2 )]} 



phy sical 



1 " £*C A {-I (l " i) 2 + (l " *) [I + | In (#)] - f + fin (=f) } 



p 2 =0 



o(g% 

(3.28a) 



15 



and 



c 2 

ghost 



dT, ghost (p) 



dp 



2 



T)2 — Yfl 

1 phy steal 



(3.28b) 

= 1- [l + I m -I + Iln (^)}| p2=Q + C% 4 )- 

We do not calculate c 2 luon because it is not used in this paper. It is important to refrain from taking 
the on-shell limit (p 2 0) until the IR divergent logarithms above have cancelled algebraically 
with similar logarithms in the amputated Green functions which are being converted to S-matrix 
elements. 

We have generalized the notion of residue to include the IR divergent terms that arise in a 
massless theory. Notice that we are extracting the residue of the propagator pole by taking the 
partial derivative of the denominator with respect to p 2 at p 2 = 0, rather than simply dividing the 
denominator by p 2 . The partial derivative extracts that piece of the IR-divergent logarithm which 
is proportional to p 2 . This piece is necessary to ensure, for example, that the two-gluon matrix 
element of the energy-momentum tensor, Eq. (4.2), is IR- finite on shell and equal to its correct 
value. 

3.7 Covariant Gluon Operator 

In [8], Hamberg and van Neerven calculate the anomalous dimension of the covariant gluon operator 
to two-loop order, C(aJ), with all free Lorentz indices contracted with a null-vector A. This selects 
the highest-spin part of the operator and eliminates the need to calculate the trace terms. 
The covariant gluon operator is 

O^-^(x) = \^-*S[F^{x)b^{x)i)i£ Ci {x) . . .D^ :i \Jx)F/™(x)] + trace terms, (3.29) 

where S denotes symmetrization of the Lorentz indices m and the trace terms make the operator 
traceless under all possible contractions of the free Lorentz indices in pairs. The c 4 - are color indices 
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in the adjoint representation. 

Hamberg and van Neerven's gauge-invariant operator is 

W = «.* A(ii>iiA(im ^ (3.30) 

where A is light-like. In units such that c = 1 = %, the mass-dimension of this operator grows 
linearly with to, 

[Cj m )] = m + 2 (3.31) 
Selecting the highest-spin piece is equivalent to selecting the lowest twist, since 

twist = mass-dimension — spin. (3.32) 

All operators of the form above, Eq. (3.30), are twist-2. 

The simplest case (to = 2) of the covariant gluon operator, Eq. (3.29), gives, up to a multiplica- 
tive factor, the gauge-invariant part of the energy-momentum tensor, Eq. (3.12) 

(n^v _ 1 fr^ip fiv _ \ n liv frpK p — _lo( G I) (r> oo\ 

a ~ 2 a pa 8^ a r pir a — 2 ' {O.OO) 

We study this case because of the relative simplicity of the calculation, but also because the 
gauge- variant operators which mix with it are supposed to be known [14,15]; they are the gauge- 
variant operators in the energy-momentum tensor, Eq. (3.13). No other operators are required. 
The gauge- variant operators for our special case also happen to be those given by the Joglekar-Lee 
prescription. We calculate the trace terms mentioned in Eq. (3.29) explicitly even though they are 
higher-twist. 

The specific case to = 2 of Hamberg and van Neerven's operator, Eq. (3.30), is 

Of = O^A.Av = \F» aP^A.A,,. (3.34) 
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Hamberg and van Neerven use a basis of operators given by Dixon and Taylor [22] before the 
BRST symmetry was fully developed. For the case m = 2, their alien operator does not correspond 
to our GV operator (Eq 3.13). The basis of operators that they chose to mix with their GI operator, 
Eq. (3.34), is 

°«L = ( P rD p ah A\ + Wd^A^ + 0(g 2 ), (3.35) 

where 



°aLn =l\FrD p ah A\ + FP"b p ab A%] + fiad^Ua + terms proportional to + O(f ) 
= h[Fr(D P A v ) a + F^iDpA^a] + fiad^Ua + terms proportional to g^ + O(f ). 



(3.36) 



This operator is not BRST-exact, in fact its BRST variation does not vanish. 

We can give a schematic form for the operators which are BRST ancestors to the alien operators 
of highest twist that mix with the GI operators for all even moments m simply by counting mass- 
dimension and twist, and by requiring an SU(N) singlet. 



an^stor (O^) = £ C H S 

2 = 



m — 2i — l 



(fj-A" 1 ) Y[(A^-A^+ 1 ) ]J (d^-k+i 

j=l k=l 



o{g), 

(3.37) 



where S denotes symmetrization of the Lorentz indices and the coefficients c\ m ^ are arbitrary (they 
are determined as a result of the renormalization). The derivatives may act on any combination of 
the fields. The dot products represent contractions over color indices. The alien operators obtained 
from these ancestors, of course, have vanishing BRST variation (modulo the trivial equation of 
motion), because of the nilpotence of that variation. 

Some examples should clarify the notation. The ancestor of the twist-2 GV part of energy- 
momentum tensor (second moment of the covariant gluon operator) in schematic form is 

an^stor = C^[fj a A^ 8^ + f, a A^d^} + O(g). (3.38) 
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Compare this with the twist-2 part of Eq. (3.14). The fields upon which the partial derivatives act 

(2) 

are not specified in the schematic form. The arbitrary constant Cq ' can be absorbed into entries 
of the renormalization mixing matrix. The O(g) terms turn out to be unnecessary. 
The dimension-6 GI operator (fourth moment of the covariant gluon operator) is 

= -U[F^b^ iC2 b^F p ^} + trace terms. (3.39) 

The BRST ancestor of the GV operators that mix with the twist-2 (highest-twist) part of the GI 
operator above is 

an^stor (OXT^) = C^S[f, a A^ d^d^d^] + C^S^A^ if if d^] + 0(g). (3.40) 

Again, the partial derivatives may act on any combination of the fields. Only one of the two 

c (4) 

arbitrary constants may be absorbed into the mixing matrix; the ratio — ^jj is determined by the 
renormalization of O s w ' i2 '' 3W . 

4 Results of the Calculation at Zero Momentum Transfer 

The results in this section were obtained under the possibly questionable operation 12 of first taking 
the limit of zero momentum transfer, before any of the Feynman diagrams are evaluated. The 
software package used was the symbolic manipulator, FORM [23], written by J. A.M. Vermaseren. 
We now evaluate Green functions of various pieces of the energy-momentum tensor first with two 
gluon fields, then with one ghost field and one anti-ghost field. 

The fields in the inserted operators are bare, while the external gluon and ghost legs are renor- 
malized as usual, since it is the renormalized fields that interact and are loop-corrected. This differs 
from Hamberg and van Neerven's calculation in which the external fields as well as the fields in the 

12 One might be suspicious of interchanging the order of the limit as the momentum transfer goes to zero with the 
other limits involved in the renormalization procedure. 
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inserted operator are bare. The two calculations therefore differ by factors of some multiplicative 
renormalization constants, Eq. (3.21). 



4.1 Green Functions of 0^ with Two Gluon Fields 



While we give the pole pieces in their entirety, the finite parts have been simplified for clarity in 
presentation. The full tensor structure can be found in the Appendix, where we also list the twist-2 
(spin-2) piece of the operators. 

4.1.1 Entire Energy-Momentum Tensor 

Consider the amputated gluon two-point Green function with the entire energy-momentum tensor, 
Eq. (3.11), inserted at zero momentum transfer. The external gluon fields have not been contracted 
with physical polarization vectors, we have not multiplied by the modified LSZ residue of the gluon 
propagator pole, and the external momenta have not been put on shell. This is what we will mean 
by an off-shell gluon Green function in the sections that follow. Explicit calculation gives 



/ 



(0|TA(j aO^i/Aj- ft \0) Amputated —PfiPi/9crT^ab 



"2 + iferCU { 



V 



ii-i 



3 + ln 



P 2 



_32 _ 10! 
9 3 111 r .2 



(4.1) 



+UV-finite terms that vanish on shell + 0(g 4 ). 



Notice that this object is UV-finite. A glance at Eq. (A. 2) in the Appendix will satisfy the reader 
that the terms not included above are also UV-finite, even off shell. The UV-finiteness supports 
the results of Freedman et al. [14, 15] on the renormalization of the symmetric energy-momentum 
tensor. Also, since there is no pole, the contribution to the anomalous dimension is zero. 

Taking into account Hamberg and van Neerven's different dimensional regularization prescrip- 
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tion (dimension 4 + e instead of 4 — 2e) and their use of bare fields instead of renormalized fields 
for the external legs, our result contracted with light-like vectors, Eq. (A. 3), agrees with theirs. 

We now do what Hamberg and van Neerven argue, quite reasonably, in their paper is impossible. 
We construct a matrix element between massless gluon states, but to do so we use the modified 
LSZ prescription described in Section 3.6. 

Contracting with physical polarization vectors, using the modified LSZ residue of the gluon 
propagator pole, c 2 gluon in Eq. (3.28a), and putting the external momenta on shell, we get the 
relatively simple S-matrix element 

{e 1 ,p,a\0 l _ tlJ \e 2 ,p,b) = -2-p tl -p l ,b ah a\ ■ e 2 + 0(g 4 ), (4.2) 

where we have used the fact that e 4 - is a physical polarization vector satisfying 

p-€i = i=l,2. (4.3) 

The physical state |e;,p, a) is meant to represent an on-shell gluon of momentum p, polarization 
vector €i, and color a. 

The modified LSZ procedure eliminates the IR divergent logarithms even before the external 
momenta are taken on shell. The result Eq. (4.2) is not surprising since 6^ is the conserved Noether 
current and 

P v = j rf 3 x Qv {x) (4.4) 

is the Noether charge. It measures the physical (non-IR- divergent) energy-momentum in a state. 
A correct calculation should show that all the higher-order corrections to the right-hand side of 
Eq. (4.2) vanish. 

In the next two sections, we examine the GI and GV pieces separately. 
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4.1.2 Gauge-Invariant Part 

Consider now the amputated off-shell gluon two-point Green function with only the gauge-invariant 
piece of the energy- momentum tensor, Eq. (3.12), inserted at zero momentum transfer. Explicit 
calculation gives 



(0|TA(j aQfiv A T ft|0) Amputated — t iq v 2 C A$ab < 



P ( 2 9<?tQiiv ~\~ 9(7h9tv ~\~ 9(71>9th) 

-\{PrPv9a^ + PrP^au + PaPv9r^ + PaP^TU ~ PaPr9^ 



t 



+PlJ,Pv9aTt>ab 



I"* 



-2 + ±,C A { 



1-* 



86 _ 10 
" 9 3 



6 + In 



In 



P 2 



+UV-finite terms that vanish on shell + 0(g 4 ). 



(4.5) 



There is a UV pole in Eq. (4.5), but this UV divergence vanishes on shell with physical polarizations. 

A further clue that we have performed the calculation correctly is the fact that Eq. (A. 4) satisfies 
the following Ward identity 



PctPt(Q\T 'A a 0,9^ A T b\®) Amputated = 0. 



(4.6) 



If we now put this result on mass-shell and use the modified LSZ procedure to derive the 
S-matrix element, we get 



I 9 2 

(e 1 ,p,a\6^\e 2 ,p,b) = P^pJab(*i ■ e 2 { -2 + ~^^ C A 



If 1 



2 ( l ~X - 31 -A I+6 



1 



f 0(g% 
(4.7) 

Notice that this physical matrix element of a gauge-invariant operator depends on the gauge-fixing 
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parameter, A. Also, the GI part is not equal to the total, so we must calculate the GV part. 



4.1.3 Gauge- Variant (Alien) Part 



Consider the amputated off-shell gluon two-point Green function with only the gauge- variant piece 
of the energy-momentum tensor, Eq. (3.13), inserted at zero momentum transfer. Explicit calcula- 
tion gives 



(Q\TA a a 0^ v A T b\0) Amputated = 7i6^2^ A^ab < 



2(1 \ 
P \29crr9liv QafiQrv QavQrfi) 

+p M tK^c a {-\ (i - i) 2 + 3 (l - i) - e} 



+UV-finite terms that vanish on shell + 0(g 4 ). 



(4.8) 



The pole terms cancel between the GI and GV parts off shell. On shell, each pole piece vanishes 
individually. 

Going on-shell and using the modified LSZ procedure, we get the S-matrix element 



9 

(ei,p,a|6^ y )|e 2 ,P,&) = P^pJab(*i ■ e 2 —r^C A 



1 



1 



3U-I 



+ 0(g*). (4.9) 



Notice that the finite part of the physical matrix element does not vanish. But, 6^7 ; is BRST- 
exact, as we observed, and there is a general theorem that BRST-exact operators have vanishing 
physical matrix elements. Thus, we know that we have a contradiction with general theorems. 

Hamberg and van Neerven have a similar result implicit in their formulae (they did not remark 
on it), but since their alien operators are not BRST-exact, physical matrix elements of their alien 
operators would not be expected to vanish. 
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4.2 Green Functions of 0^ with One Ghost Field and One Anti-Ghost Field 
4.2.1 Entire Energy-Momentum Tensor 

We have just seen that the two-gluon matrix element of the gauge- variant part of the energy- 
momentum tensor is nonzero. Since 6$ V ^ is BRST-exact, this contradicts a crucial part of the 
theory on the renormalization of gauge-invariant operators and so we cannot take for granted any 
of the results of this theory, but must verify the results. 

In particular, we need to verify the finiteness at one-loop order of Green functions of the energy- 
momentum tensor. The previous section has established this for the gluon two-point Green function 
and in this section we verify finiteness for the ghost-anti-ghost Green function. 

All the necessary counterterms are determined by the formula for 0^; they are obtained by 
expanding 6^ in terms of renormalized fields by Eq. (3.16). 

Consider the amputated off-shell ghost two-point Green function with the entire energy-momentum 
tensor, Eq. (3.11), inserted at zero momentum transfer. Explicit calculation gives 



(O\Tu) a lll/ 'l] b \O) Amputated ab(2p 'fiPv ~ P 2 9 ^ 



A ) 1 2 



-1 + In 



1 +ln 



P 2 



P»Pv - 4 m ( — ? V 



1-iln 



+0(g 4 ) 



(4.10) 



Using the modified LSZ residue of the ghost propagator pole, c 2 ghost in Eq. (3.28b), and putting 
the external momenta on shell, we get the comparatively simple S-matrix element 



(p, a 1 0fj, u \p,b) = 2p^pJ ah + (D(g 4 ), 



(4.11) 



which is correct for the expectation value of in a properly normalized state of momentum p. 
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Here the state vector \p, a) is meant to represent an on-shell ghost of momentum p and color a. 
Again, although we only performed the one-loop calculation, all higher-order corrections should 
vanish. 

The twist-2 (spin-2) piece of the amputated Green function Eq. (4.10) above in which the free 
Lorentz indices of the inserted operator, fj, and is, are contracted with a null- vector, A, is 



(4.12) 



(OlTLJaA^ 9 ^ A u ?76 1 0) Amputated = 

(p-A) 2 6 ab (2+^C A ^(l-{) l+ln[^) -l + ln 
4.2.2 Gauge-Invariant Part 

Consider now the amputated off-shell ghost two-point Green function with only the gauge-invariant 
piece of the energy-momentum tensor, Eq. (3.12), inserted at zero momentum transfer. Explicit 
calculation gives 



(0 1 Tu a 0^ Vb\0) Amputated =\-^C A b ah {p ll p v - \p 2 CJ ^) 



16-7T 



-CA^ab 



1 - In 



P 2 



Py.Pv 



P 2 ) 2 

— P 9^ 



o(g 4 ). 

(4.13) 



The twist-2 (spin-2) piece of this amputated Green function in which the free Lorentz indices 
of the inserted operator, fj, and is, are contracted with a null- vector, A, is 



(0\Tu a A^l G ^A" Vb \0) Amputated = (p ■ A) 2 -^-C A 6 ab 



i + l-ln (=4 



(4.14) 



4.2.3 Gauge- Variant (Alien) Part 

Consider the amputated off-shell ghost two-point Green function with only the gauge- variant piece 
of the energy-momentum tensor, Eq. (3.13), inserted at zero momentum transfer. Explicit calcula- 
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tion gives 



(0|Tw a 6»;L 'Vb 1 0) Amputated = 7if^ C ' A& ab{~P ftPv + \p' '9 >) 



1-i 1 



1 +ln 



+2 



-1 + In 



Py.Pv 



PlcPk ~ \ In 



l-|ln 



A 2 

— P 9^ 



P 2 9^ 



+o(g 4 )- 



(4.15) 



Notice that the pole terms cancel between the GI and GV pieces. 

The twist-2 (spin-2) piece of this amputated Green function in which the free Lorentz indices 
of the inserted operator, fj, and is, are contracted with a null- vector, A, is 



(o|r^ a A^^ y) A^ 6 |o) Amputate(i = 



1 + In 



2 + 2 In 



0(9 A ) 
(4.16) 



4.3 Renormalization Mixing Matrix 

If we do not require a priori that the matrix be triangular, then the most general form is 



r \e 



R 



R [Efj 



Zgg Zqa Zqe 



Zag Zaa Zae 



Zee 



\ ( \ 

} (GI) 



q{GV) 



(4.17) 
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where the operator of class E, which vanishes by the equations of motion (Eqs. 3.6) and mixes with 
the operators in the energy-momentum tensor, is 

=S ab A^ c [(b p FP») c + \d v d ■ A c + gc ade (d»f, d )u e ] 
+S ab A v c[(b P F^) c + \d»d ■ A c + gc cde (d^f, d )Cj e ] 

-\Kb9^K c[(b p FP*) c + Xd^d ■ A c + gc cde (d*r, d )u e ] ( 4 - 18 ) 

+ \£>abg^[(b p d p f}) c ]Lb c a 

+ l^abg^flc[d p (b p cb) c ] (1 - a). 

Like the fields in the energy-momentum tensor, the fields in the operator E^ v above are bare. The 
parameter a in the last two lines above is free to vary between and 1 since the matrix elements 
considered are not sufficient to distinguish between the equations of motion for the ghost and 
antighost fields. 

We find that the following elements of the mixing matrix are compatible with both the two-gluon 
and two-ghost projections: 

Z GG = 1 + 0(g 3 ) 
ZGA=-\\^C A + 0(g3) 
ZGE=\\^C A + 0(g3) 

Zag= 0(g 3 ) ( 4 - 19 ) 
ZAA=l + \\^C A + 0(g3) 
ZAE=-\\^C A + 0(g3) 
Z EE = 1 + 0(g 2 ). 
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The renormalization mixing matrix is triangular to 0(g 2 ). We do not calculate Zee beyond the 
tree level. 

4.4 BRST Ward Identity (Slavnov- Taylor Identity) 

We have seen that the gauge- variant part of the energy-momentum tensor is nonzero in an on-shell 
matrix element. However, the gauge- variant part of the energy-momentum tensor is BRST-exact, 
in accordance with general theory and a very simple proof states that physical matrix elements of 
such operators vanish [13]. 

In this section, we resolve the contradiction. The proof that physical matrix elements of BRST- 
exact operators vanish proceeds by using a Ward identity based on the BRST variation to relate 
the matrix element under study to a particular Green function of the ancestor operator of the 
BRST-exact operator. This Green function has a manifest factor of zero when put on shell, but 
the zero is compensated by a quadratic infra-red divergence present only when the matrix element 
of the ancestor operator is evaluated at zero momentum transfer, as we will now see. 

Let us assume that the BRST Ward identity holds for unamputated Green functions off shell, 
calculated at zero momentum transfer. 

The BRST variations of the bare fields are given in Eq. (3.7). We need the BRST variation of 
the renormalized fields in terms of renormalized fields and parameters. This is sometimes referred 
to as the 'renormalized BRST variation' and is related to the canonical BRST variation by factors of 
the renormalization constants, Eq. (3.21). The goal in defining a renormalized BRST variation is to 
obtain UV finite Green functions with renormalized fields. The renormalized constant Grassmann 
parameter which accomplishes this goal [13] is 



i i 



8Z = Z A *Zr,HZ 
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so we have 



z% ) ^ 



- s r 

<WsT^a = -^ZgZigCabcUbUJC =>- B ^ bT ^a= ~ 2 %g Z Z\gC a b c LObLO c 



(4.21) 



A n - ZxZ A \fi . a T? ^BRST — \f)' A 

°BRST'/a- 1 A(J SiaVQ r~ l/a— AU A a 

s r s r s r 

The three operators - B f^ L A l _ t a , - B f^ L uj a , and - B f^- L i] a are all finite. 

The finite ancestor operator is 

ancestor = (d^A^ a + {d^A, a - g llv [\r la d ■ A a + (d p rj a )A p a ] , (4.22) 

where 

^ancestor = 9^V) _ Ig^d^D^ . (4.23) 

Remember that the fields in the energy-momentum tensor operator are bare. The bare ancestor 
and the finite (renormalized) ancestor are related by 

ancestor (o^j v ^ = Z\z\ ancestor (flg 1 ^) (4.24) 

Now, the BRST variation of any Green function vanishes. Consider the particular case 

S BRST (0\TA a ancestor (0<f/ ) ) A T 6 |0) = 0. (4.25) 
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This gives 



= (0 



T[<5 BR stAt a]ancestor (o^ V ^ A T b 



TA aa \b BBST ancestor {o^ V ^ A rb 



TA a a ancestor (o^^) [<$ B rst^t b] 



(4.26) 



The vanishing of Eq. (4.25) and the chain rule for the BRST variation can be proven by defining 
the variation in terms of (anti)commutators with the Noether charge associated with the BRST 
symmetry [19]. 

Equation (4.26) becomes 



=( 



T 



ZodtjLUa - Z g Z Q Z\gc ade A a d u e S£ 



ancestor ( 6^ ' ) A T b 



o\TA aa [0^ V) S^\ A Tb \o 



TA„ a \-y^fj c dP(D p u) c 6S\ A Tb 



TA<j a ancestor [ 6$ V ^ 



Z d T uj b - Z g Z Z 2 A gc bde A T d cj e )S^ 



(4.27) 



The constant Grassmann parameter, 8£, can be factored out of Eq. (4.27) above if it is anti- 
commuted through the ancestor operator which has a Grassmann parity of 1 (because each term 
contains one antighost field). This is responsible for the relative minus sign below. 
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T 



Zod a Lo a - Z g Z Z\gc ade A a d L0 e 



( GV ) 

ancestor ( d\ v ' ) A T b 



TA a ji G / ] A Tb 



TA aa \ -\g^r) c d p (b p Cj) c \ A Tb 



TA<j a ancestor ( 6$ V ^ 



Z d T oj b - Z g Z Z\gc bde A T d L0 e 



(4.28) 
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One must also keep in mind that the order of the Grassmann fields, f) (in the ancestor operator) 
and u, in the last two lines above is opposite to the canonical ordering and that 



{0\Tf)Uj(Operator)\0) = — {0\Tujf)(Operator)\0) 



(4.29) 



It is now obvious that there exists an alternate calculation which will provide the two-gluon 
physical matrix element of 9^ ' . 



T Z d a oj a - Z g Z ZXgc ade A a d u e 



TA aa \-\g^fi c d p (b p Lo) c \ A Tb 



ancestor 



TA<j a ancestor [9^ V ^ 



Z d T oj b - Z g Z ZXgc bde A T d u, 



(4.30) 



The usual proof that this matrix element vanishes on shell relies on the assumption that the 
triangle graph depicted on the bottom of page 46 in the Appendix does not contain a quadratic 
infra-red singularity. This graph and another like it contain an unusual vertex, one gluon and 
one ghost at the same space-time point. There is no external line, therefore nothing to amputate, 
but the amputation procedure for other graphs, which have ghost legs, requires that all diagrams 
be multiplied by the inverse of a ghost propagator, which of course is proportional to p 2 . If the 
diagrams contain at worst logarithmic IR divergences, then the additional factor of p 2 will cause 
them to vanish on shell. The diagrams with an external ghost line have a derivative acting on 
the ghost field. In momentum space, the derivative becomes a factor of p a which gives zero when 
contracted with the physical polarization vector, e a , associated with the gluon leg. 

Explicit calculation of the relevant graph, with the result Eq. (A. 9), shows that there is, in 
fact, a -p- divergence which cancels against the inverse of the ghost propagator introduced in the 
amputation procedure, as we claimed at the beginning of this section. This IR pole occurs only at 
zero momentum transfer. 



31 



5 Conclusion 



We have seen by explicit calculation that one of the central results of Joglekar and Lee, that physical 
matrix elements of BRST-exact operators must vanish, fails at one loop order. We give a form for 
the BRST ancestors of the alien operators required in the renormalization of the covariant gluon 
operator. Our alien operators are then manifestly BRST-exact whereas the basis of alien operators 
proposed by Dixon and Taylor (those used by Hamberg and van Neerven) are not BRST-exact. The 
Dixon and Taylor set of alien operators are not guaranteed to vanish in physical matrix elements. 

We have verified the predictions of Freedman et al. on the finiteness of the energy-momentum 
tensor in both gluon and ghost two-point functions to one loop order at zero momentum transfer 
by evaluating diagrams with a BRST-exact alien operator insertion. 

The BRST Ward identity demonstrates where the proof of the Joglekar-Lee theorem breaks 
down. Taking the momentum transfer to zero too soon introduces spurious infra-red divergences 
which cancel factors of zero on which the proof relies. Calculations performed using the Dixon and 
Taylor set of alien operators cannot by analyzed through the BRST Ward identity. 

The physical region of interest in almost all calculations involving the renormalization of compos- 
ite operators, such as the calculations required in the operator product expansion, is the exceptional 
point of zero momentum transfer. To expedite the computation, one sets the momentum transfer 
to zero at the very beginning, thereby eliminating one scale from the problem. In some calculations 
involving final state cuts, it is not clear how one would generalize to nonzero momentum transfer. 

The alternative is to keep the momentum transfer arbitrary until after the Feynman graphs 
have been evaluated, and only then to set the momentum transfer to zero. With this procedure, 
the Joglekar-Lee theorem should apply, making it unnecessary to compute the graphs containing 
the alien operator insertion. The price to be paid, of course, is the introduction of another mo- 
mentum scale into the problem and a corresponding increase in the complexity and volume of the 
analysis. However, it is not obvious that the limit of zero momentum transfer is the only source 
of contradiction with theory. A recent calculation by Harris and Smith [16] at nonzero momentum 
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transfer indicates that the discrepancy persists. 

There has been a sense of disquiet in the literature about zero momentum transfer for a long 
time. Joglekar mentions in the concluding section of [24] that, at the exceptional momentum point, 
Q = 0, matrix elements of gauge-invariant operators lose some of the properties that make them 
manageable at nonzero momentum transfer. C. Lee [25] works with the twist-2 piece of the energy- 
momentum tensor to show that certain pieces of the calculation at zero momentum transfer can 
yield useful information, that is the coefficients of certain terms are the same, independent of the 
momentum transfer. He calculates only the pole terms of the two-gluon Green function at one-loop 
order at both zero and nonzero momentum transfer. The unease was certainly justified; some results 
hold while others fail utterly. It is not unreasonable to question all calculations performed when 
the limit of zero momentum transfer was applied initially, and such calculations are the mainstay 
of perturbative QCD. 
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Appendix 

A.l Right Derivatives 

Right derivatives [26] are such that 




(A.l) 
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where P Y is the 'Grassmann parity' of the quantity Y . The (anti)ghost field components have 
Grassmann parity 1, while the c-number parameters and boson field components have Grassmann 
parity 0. Fermion field components would be assigned Grassmann parity 1. 

A. 2 Complete Off-Shell Calculations for Two-Gluon Amputated Green Func- 
tions 

In this part of the Appendix, we give the full Lorentz structure for the amputated Green functions 
of the energy-momentum tensor operators with two gluon fields off mass-shell, at zero momentum 
transfer, to one loop order. The ^ poles are purely UV divergences while all of the IR divergences 
(as p 2 0) are seen as logarithms. 
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A. 2.1 Entire Energy-Momentum Tensor 

{0\TA a a^iiuA-r b\®) Amputated = <*> ablp' \ ( JutQ 'fiu ~ 9<J{i9tv ~ 9<jv9tv) 
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(A.2) 
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To make the comparison with Hamberg and van Neerven more transparent, we present the 
twist-2 (spin-2) piece of this amputated Green function in which the free Lorentz indices of the 
inserted operator, ji and z/, are contracted with a null-vector, A: 
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(A.3) 
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A. 2. 2 Gauge-Invariant Part 
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The twist-2 (spin-2) piece of this amputated Green function in which the free Lorentz indices 
of the inserted operator, fj, and is, are contracted with a null- vector, A, is 
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(A.5) 
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A. 2. 3 Gauge- Variant (Alien) Part 



{Q\TA a aO^" 1 A T b\0) Amputated = ~ 7 -^iCA^ab [p 2 ( ~ \ ( JaT ( Jiiv + 9aii9ru + 9<jv9t 
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p 2 (g 
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(A.6) 
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The twist-2 (spin-2) piece of this amputated Green function in which the free Lorentz indices 
of the inserted operator, fj, and is, are contracted with a null- vector, A, is 



{0\TA a a A"0g V) A"A T 6 |0) AmV utated = ~\^C A S ab [2p 2 A a A T - (p ■ A)( Pa A T + p T A a ) 



-2\S ab (p-A)( Pa A T + p T A a ) 
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(A.7) 
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A. 3 Unamputated Two-Gluon Green Function of Alien Piece of the Energy- 
Momentum Tensor Off-Shell 

In order to use the BRST Ward identity (which is valid off mass-shell), we need the unamputated 
Green function which we obtain from the amputated Green function, Eq. (A. 6), by attaching the 
dressed external gluon propagators, Eq. (3.24a). 13 

(O\TA^jf/ ) A Tb \O) = V^ l '( P )(O\TA a , a ,0f/ ) A T , bl \0) Amp utatedV^(p) 



1-* 
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+o(g 4 ) 



(A.8) 



13 The modified LSZ prescription, rather than the dressed external propagators, would have been used to go on 
shell. 
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A. 4 Feynman Diagrams 

The following graphs contribute to the matrix elements at one-loop order. The external propaga- 
tors are amputated; one-loop corrections to the legs are handled by the modified LSZ reduction 
procedure explained in Section 3.6 to go on shell or by attaching dressed external propagators as 
in Section A. 3 to remain off shell. The inserted composite operator is represented by the symbol 
(g). In this section, (g) = 6^, offi\ or 6^ V \ 

A. 4.1 Gluon Two-point Function 



The Born graph: 




The order g 2 graphs: 
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The following diagram vanishes for zero momentum transfer, Q, but contributes in the nonzero 
momentum transfer case: 



it <?#o 




The following diagram always vanishes in dimensional regularization, regardless of the momen- 
tum transfer, because the integral contains no momentum scale: 




VQ 
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A. 4. 2 Ghost Two-point Function 

The Born graph: 




The order g 2 graphs: 




A.4.3 The BRST Ward Identity Graphs 

In this section, the symbol ® stands for the renormalized BRST ancestor of 6^ , (Eq. 4.22). 
The Born graph: 




The order g 2 graphs: 




The two graphs above have a derivative acting on the ghost field at the end of the external line. 




The last two graphs contain an unusual vertex (a vertex not contained in the Lagrangian 
density), one ghost and one gluon field at the same space-time point. Also, they are each multiplied 
by the inverse of a ghost propagator as part of the amputation procedure. The three graphs above 
vanish in physical matrix elements: the two graphs with external derivatives will be proportional 
to the external momentum and this contracted with the physical gluon polarization will give zero; 
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the graph immediately above contains a free Lorentz index in its unusual vertex and must be 
proportional to the external momentum as well. Only the graph below contributes to on-shell 
matrix elements. 
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(A.9) 



{^9ct)i9tv t 9cw9t)i) 



n PnPv _3 Pffji PvPn - 

9 ff(TT 9 1 9 S'r/i T 9 SVi/ 9 9 

P \^ 



+ (?G/ 4 ) 



This graph plus its mirror image (a <r+ r, a <r+ b, p <r+ —p) contribute to the two amputated 
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Green functions 



T 



TA<j a ancestor [ 6$ V ^ 



ancestor ( 6^ ' ) A T b 



' Amputated 



(A.10) 



Amputated 



corresponding to the unamputated Green functions found in Eq. 4.30. The terms proportional to 



V 



—g aT in the finite part ruin the proof that the physical matrix element of an alien operator must 



vanish. 14 

To isolate the parts that survive on shell, we contract Eq. (A. 9) with a physical gluon polariza- 
tion vector, e T , to obtain 
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(A.11) 



+3-^e„ - £ ( — —e„ H ^e, 



14 The pole piece cannot contain a -p divergence because of locality. 
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A. 5 Feynman Rules 

In this section, we give our conventions for the Feynman rules of common objects in pQCD and 
the non-standard vertices introduced in this paper. 

A. 5.1 Propagators 



/i,a z/,b 



te 



9 iip 



1 - - 



1 \ k,,k 



XJ k 2 + ie\ 



(A.12) 



k 2 + ie 



(A.13) 



Include a factor of —1 for every ghost loop. 



A. 5. 2 Vertices in the Lagrangian Density 



All momenta are defined to flow into the vertex under consideration. 




-gc aia2a3 [g^( Pl -p 2 p + g™(p 2 - P3 p +g^(P3 - PiP] (A.14) 



z/i,ai 




z/4, 04 



-i g^cba^cba^ig^g^ - g^g^) 



+ Cba 1 a 3 Cba i a 2 (9 l/ll/i g l/3l/2 ~ </^</^) 
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(A.15) 
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/i,a 



A. 5. 3 Non-Standard Vertices at Zero Momentum Transfer 



(A.16) 
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E 



(A.22) 
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